Field Driven Pairing State Phase Transition 
in d x 2_ y 2 + zda^-Wave Superconductors 
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Within the framework of the Ginzburg-Landau theory for d x 2_ y 2 + id xy -w&ve superconductors, 
we discuss the pairing state phase transition in the absence of the Zeeman coupling between the 
Cooper pair orbital angular momentum and the magnetic field. We find that above a temperature 
T t , the pairing state in a magnetic field is pure d x 2_ y 2-wa,ve. However, below T», the pairing state 
is d x 2_ y 2 + idajy-wave at low fields, and it becomes pure d x 2_ y 2-wicve at higher fields. Between 
these pairing states there exists a field driven phase transition . The transition field increases with 
decreasing temperature. In the field-temperature phase diagram, the phase transition line is obtained 
theoretically by a combined use of a variational method and the Virial theorem. The analytical 
result is found to be in good agreement with numerical simulation results of the Gingzburg-Landau 
equations. The validity of the variational method is discussed. The difference to the case with the 
Zeeman coupling is discussed, which may be utilized to the detection of the Zeeman coupling. 
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The phase-sensitive experiment with a tricrystal super- 
conducting ring magnetometry jjj demonstrated that in 
high temperature superconductors, the dominant pair- 
ing channel is the d x 2 _ y 2 -wave one. However, it is not 
yet clear whether there were some sub-dominant pairing 
channels. It is even less clear what would be the symme- 
try of a sub-dominant pairing channel if it existed at all. 
Such questions arise from a number of experiments, e.g., 
the observation of surface-induced broken time-reversal- 
symmetry (T hereafter) in YBCO tunnel junctions [Q, 
the observation of fractional vortices trapped in a bound- 
ary junction ||, and the abnormal field dependence of 
the low temperature thermal conductivity n e in BSCCO 
superconductors [QJ^]. These unusual phenomena can 
not be adequately explained by the d x 2_ y 2-wetve pair- 
ing channel alone. Some sub-dominant channels (such as 
s— or d xl) -wave ones) might have played a role in these 
phenomena. Thus it is interesting to study theoretically 
the properties of the superconductors in the presence of 
sub-dominant pairing channels. 

As a model study, we consider the relevant singlet sub- 
dominant pairing channel to be the d xy - channel, which 
has been hotly discussed recently in the context of the ab- 
normal thermal conductivities in BSCCO superconduc- 
tors. {§ [?]]. In this paper, using a Ginzburg-Landau (GL) 
theory for d x 2_ y 2 + id xy -w&ve superconductors, jTj we 
discuss the pairing state phase transition driven by the 
magnetic field. In the absence of the Zeeman coupling 
between the Cooper pair orbital angular momentum and 
the magnetic field (see below), we find that above a 
temperature T„, the pairing state in a magnetic field is 
pure d a: 2_ y 2-wave. However, below T*, the pairing state 
is d x 2_ y 2 + id xy -w&ve at low fields, and it becomes pure 
c? 2 .2_ 1/ 2-wave at higher fields. There exists a field driven 
phase transition between these pairing states. The tran- 
sition field increases with decreasing temperature. In the 



field-temperature phase diagram, we are able to obtain 
the phase transition line using a variational method and 
the Virial theorem. The analytical result is in good agree- 
ment with that obtained from numerical simulation of the 
Ginzburg-Landau equations. 

The GL free energy of a d x 2_ y 2 + id xy -wave super- 
conductor can be obtained from the modified Bardeen- 
Cooper-Shricffcr (BCS) gap equation, the Gor'kov theory 
or path integral formulation (in the weak coupling limit), 
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T[3\D\ 4 / 

+ \D\ 2 \D'\ 2 /4+{D*D' + c.c.) 2 /8] 

+K [\UD\ 2 + \UD'\ 2 )] + [ (V X A) 2 /8tt, (1) 

Jn 

where J n denotes integration over the afe-plane, IT = 
— iV — 2eA/hc is the gauge invariant gradient and 
B = V x A is the local induction (in the z-, or c-, direc- 
tion). Here D and D' are order parameters in the d x 2_ y 2- 
and djy-channels, respectively, on = N(0)lnT/Ti (i — 
D, D'), where N(0) is the normal state density of states, 
and Ti is the bare superconducting critical temperature 
in the i channel. K and T are material-dependent pa- 
rameters derivable from microscopic theories and are 
assumed to be temperature independent here for sim- 
plicity. We assume that Td> <C Td, i.e., the dominant 
pairing channel is the d;j,2_ y 2-wave channel, as this would 
be consistent with high temperature superconductors. In 
this paper, we do not consider the Zeeman coupling term 
F z cx -iB{D*D' - c.c.) in the free energy. |7|-[| On one 
hand, F z needs further identification by a strong cou- 
pling microscopic theory. On the other hand, the pairing 
states in a magnetic field in the absence of F z is interest- 
ing in its own right and should be compared to that in the 
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presence of F z , which will be discussed in the concluding 
section. 

In order to see the relative importance of the various 
contributions to the free energy and for later convenience, 
we rewrite the free energy in terms of dimensionless quan- 
tities as, 

F = E c J{-\d\ 2 - a\d'\ 2 + (|d| 4 + |d'| 4 )/2 + |d| 2 |d'| 2 /3 

+ (d*d' + c.c) 2 /6 + \nd\ 2 + \nd'\ 2 + K 2 b 2 }, (2) 

where E c = H 2 ^ 2 /4tt with H c and £ being the ther- 
modynamic critical field and coherence length, respec- 
tively, in the absence of the d xy -w&ve channel. All quan- 
tities under the integration symbol are now dimension- 
less: a = a D> /a D = \n{T/T D >)J\n{T/T D ), d = D/D , 
d' = D' /Dp, r = R /(j, a = A/A , and b = B/B . Here 
Dq = y/—4ct£)/3r is the value of D in the absence of the 
order parameter D and the magnetic field, k is the bare 
GL parameter in the absence of D', and Aq = <I>o/27r£ 
where $o is the flux quantum. Bo = Aq/£ is the cor- 
responding upper critical field. Finally, n = — iV — a 
now denotes dimensionless gauge invariant gradient (i.e., 

tt = no. 

The important difference between the above GL free 
energy Eq.(|2|) and that for a d + is-wave superconduc- 
tor lies in the fact that in the former, we have no mixed 
gradient coupling between the two order parameters. JTcj ] 
Therefore, we can expect that the single vortex remains 
to be circularly symmetric. This should be contrasted to 
the four- fold and two-fold symmetric vortices in d + is 
wave superconductors. |ll| l An even more subtle conse- 
quence of the absence of the mixed gradient is that there 
are three possible solutions, namely: the pure d solution, 
the pure d! solution, and the mixed-wave solution. The 
physical solution of the system corresponds to that with 
the lowest free energy. In this respect, it is not impossible 
that there would be a phase transition between these so- 
lutions in the field-temperature phase diagram. Indeed, 
this picture was conjectured recently from the behaviors 
of the superconductor at zero field and the upper critical 
field, § which we summarize for completeness as follows, 
(i) At zero field, it is easy to see that the two order pa- 
rameters develop a relative phase difference of ±7r/2 in 
order to minimize the free energy. This d x 2_ y 2 ± id xy - 
wave pairing state is T-breaking, which was recently ar- 
gued to be relevant to the abnormal thermal conductiv- 
ity in BSCCO superconductors. The amplitudes of 
the two order parameters are |<fp = 3(3 — oe)/8 and 
\d'\ 2 = 3(3a — l)/8. Thus, d^y-wave appears only when 
a > 1/3, with a T-breaking transition point at a* = 1/3. 
According to the above definition of a, this amounts to 

a transition temperature T* = \l T' D ^ /Tp. The zero-field 
T-breaking transition is a second-order phase transition, 
in that d' emerges continuously at T < T». (ii) From 



the linearized GL equations, it was found that the up- 
per critical field b C 2 = B c ijBo = 1 as long as a < 1 (or 
T > OK). Moreover, at the upper critical field the eigen 
solution to the linearized GL equations indicates that the 
order parameter d! vanishes identically, although it would 
be nonzero at the same temperature (< T*) but at zero 
field, (iii) Combining both aspects we believe that there 
should be a field-driven pairing-state phase transition at 
T < T* (or a > a*), where at zero and low fields, the 
pairing state is the d x 2_ y 2 + id xy -w&ve, while at higher 
fields, it transforms to the pure d x 2_ y 2-w&ve. The mech- 
anism is clear. Since the two order parameters are not 
coupled by mixed gradient terms, they are frustrated in 
the vortex state: while the coexistence of them lowers the 
homogeneous energy at a > a*, winding of both order 
parameters due to vortices increases the kinetic energy. 
The competing energies should drive a phase transition. 
Moreover, the transition field should increase with a as 
the gaining of homogeneous energy increases relatively. 
On the other hand, at T > T* the pairing state is always 
the pure d x 2_ v 2-w&ve at any fields below the upper crit- 
ical field. This picture will be discussed in detail, in the 
proceeding section. 

In our case, there are two order parameters. Even ap- 
proximate solutions to the GL equations, such as for the 
conventional superconductors, jl2| are difficult to obtain. 
Thus we shall restrict ourselves to the variational treat- 
ment of the system. The basic idea is as follows. We 
employ a reference conventional system with one order 
parameter tp, and set d and d! to vary in space in the 
same manner as ip does, but with different amplitudes 
which are to be optimized. Consequently, the behav- 
ior of the present system is then approximately given by 
these amplitudes. 

For our purpose, let us recapitulate some known results 
for the reference system. The free energy is written as 
F = E c f r {-\ip\ 2 + \tp\' L /2 + \Trip\ 2 + K 2 b 2 }. Two identities 
follow immediately from the solution that minimizes this 
free energy. Firstly, 

e 2 + 2e 4 + e fc = 0, (3) 

where e 2 = (-|^| 2 ), e 4 = (|V>| 4 /2) and e k = (\n^\ 2 ). 
Here (•) denotes spatial average. Second, 2n 2 h{b) — + 
2k 2 (b 2 ), which is known as the Virial theorem. [ p^|JT^ ]In 
the limit /« > 1, the magnetic induction can be safely 
treated as being uniform at h h c \ (with h c \ being the 
dimensionless lower critical field) , so that an approximate 
Virial theorem can be written as 

e k w 2n 2 [h(b) - b]b, (4) 

where h(b) is the magnetization curve which can be read 
off from the textbook. 113] 

We now come back to the d x 2_ y 2 + id xy -w&ve super- 
conductors. As mentioned above, we set d = jiifj and 
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d! = ±iisip, with two real and positive variational am- 
plitudes fi and v. The d! order parameter adopts a 
residual relative phase to d in order to lower the free en- 
ergy Eq.(^). Substitution of these ingredients into Eq.(2) 
yields a free energy in terms of jj,, v, and (i — 2, 4, k). 
In dimensionless form, the free energy density is 

/ = A* 2e 2 + av 2 e 2 + (m 4 + ^ 4 )e 4 

+2/xVe 4 /3 + (/i 2 + v 2 )e k + n 2 b 2 . (5) 

Minimizing / with respect to fi and v, we get: 

/i 2 = 3[(a-3)e 2 -2e fe ]/(16e 4 ); (6) 
v 2 = 3[(1 - 3a)e 2 - 2e fe ]/(16e 4 ). (7) 

At low and intermediate fields, the magnetization 
curve in the reference system in the limit n 3> 1 is: h = 
& + ln(l/6)/2K 2 (in dimensionless form). [[L2| In combina- 
tion with Eq.(Q), we have e k = 61n(l/6). On the other 
hand, in the present field regime, the distance between 
two vortices is much larger than the coherence length, 
hence we can assume e 2 — — 1 and e 4 ~ 1/2. This approx- 
imation would certainly violate Eq. (|J) up to the order of 
efc (<C 1), but suffices to yield an order of magnitude esti- 
mation of fi and v. Substituting the approximate e^'s into 
Eqs.([§ and @), we get \i 2 « 3(3/2 - a/2 - 61nl/6)/4, 
and v 2 k. 3(-l/2 + 3a/2 - 6 In 1/6) /4. We clearly see 
that /j, 2 > as long as a < 1 (i.e., d xy -chaime\ is sub- 
dominant). This means that the d x 2_ y 2-w&ve order pa- 
rameter d is always present in the system. The d xy -w&ve 
order parameter d' is present provided that v 2 > 0. Thus 
the field driven phase transition is determined by v 2 = 0, 
or 

& In 1/6= (3a-l)/2. (8) 

At high fields closer to the upper critical field, the vor- 
tices are densely distributed, hence the system can not be 
treated as above because the order parameters are dras- 
tically suppressed by the magnetic field. Fortunately, in 
this case we can resort to the high-field Abrikosov vortex 
lattice solution. In the reference system one has 2e 4 = 
f3 A e 2 , where (3 A is the Abrikosov constant. We combine 
this with Eq. (|) to find e 2 = (^1 - 4(j A e k -l)/ (2(3 A ). On 
the other hand, the dimensionless magnetization curve 
in high fields is h = b + (1 - b)/(2n 2 f3 A ), @ substi- 
tution of which into Eq.(|J) gives e k = b(l — b)/j3 A . 
Putting the above together, we have in the high field 
regime: fi 2 = 3[(a — 3)e 2 — 2e k ]/ (Sf3 A e\), and v 2 = 
3[(1 — 3a)e 2 — 2e/ ! ]/(8/3 J 4e 2 ! ), where both e 2 and e k can be 
expressed as functions of b. Again the condition v 2 = 
determines the field-driven phase transition line: 

b = (1 + V9a 2 - 12a + 4)/2. (9) 

The phase transition lines Eqs.(||) and (||) are the main 
results of this work. They are plotted in Fig.l (solid 



lines) . Note that we have skipped the unphysical portions 
of Eq.(||) at high fields and Eq.@ at low fields. Also note 
that the dimensionless upper critical field is unity, but the 
physical upper critical field Bo depends on temperature. 
On the other hand, a depends on temperature also (see 
the definition above). The field regime between the two 
solid lines in Fig. 1 can not be determined theoretically by 
our approach, in that the order parameters are already 
suppressed by the field so that the London approximation 
fails on one hand, but the Abrikosov solution is still not 
reliable enough on the other hand. 

In order to check the accuracy of the variational 
method, we now obtain the transition line by numeri- 
cal simulation of the GL equations derived from Eq.(2). 
The simulation is performed in a unit cell of the vortex 
lattice and the magnetic induction is treated as being 
uniform for simplicity. The latter assumption is suitable 
at not too low fields and at n 1, which is relevant to 
high-T c superconductors. The simulation method is well 
documented in the literature. [l5| From our simulation 
result, it is verified that the local relative phase differ- 
ence between the two order parameters is indeed ±7r/2. 
This result provides a strong support to the validity of 
our analytical treatment. The maximum amplitude of 
the order parameter d' in the vortex solution, namely, 
I I max = | -D' | max /-Do j as a function of a at various mag- 
netic fields is plotted in Fig. 2, where we can clearly see 
that at a fixed field, d 1 drops to zero at a specific value of 
a. This signals a field-driven phase transition from the 
d x 2_ y 2 + id xy - wave pairing state to the pure d^2_ y 2-wave 
one (with vortices in the system). The transition value 
of a decreases with decreasing magnetic field. The set of 
transition points are plotted in Fig.l (squares). In the 
high field regime, the analytical result is in good agree- 
ment with the numerical result, whereas in the low field 
regime, the result is only in qualitative agreement. This 
is understandable from the fact that we have neglected 
the vortex core energy and have adopted a crude approx- 
imation for the spatial variation of the order parameters 
in our analytical treatment. In principle, the two order 
parameters have different coherence lengths. They de- 
termine the length scales for the spatial variations of the 
order parameters. Our approximation is equivalent to as- 
sign the same coherence length to both order parameters. 
The good agreement at high fields is also understandable 
from the view that in this case the vortices are strongly 
overlapped so that the length scale approximation is not 
essential. In numerical simulations we find some signs of 
abrupt drop of d' as a function of b or a in the parameter 
space. This may points to a weakly first-order transition. 
However, our analytical result is clearly a second order 
phase transition within the specified approximation. 

Before closing, let us comment on the difference of our 
system to that with a Zeeman coupling. In the latter 
case, it was shown that a d x 2 _ y 2 + id xy pairing state is 
always induced by vortices. In other words, there 
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is no further pairing state phase transition in interme- [17] Mei-Rong Li, P. J. Hirshfeld, and P. Wofle, cond 



diat p magnetic fields. Moreover, Balatsky [ 16 1 recently 



argued that the Zeeman coupling enhances superconduc- 
tivity near the upper critical field. These distinct differ- 
ences to the case discussed above can be utilized for the 
detection of the Zeeman coupling. 
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FIG. 1. Solid lines: analytical phase transition line; 
Squares: transition points extracted from numerical simula- 
tions. The error bars indicate the increment of the field in 
the field-scanning. 
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FIG. 2. The maximum of \d'\, i.e., -D'|max/-Do as a func- 
tion of a at various fields. The arrows indicate the estimated 
positions at which d! drops to zero. 
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